scUnifrac is to quantify cell subpopulation diversity between two single-cell transcriptome profiles, which calculates the distance, estimates the statistical significance of the distance, identifies the subpopulations that are significant different between two profiles, finds genes that mark subpopulations, and predicts cell types of subpopulations. If two single-cell RNA-seq profiles have identical cell subpopulation structures (the same subpopulations with the same proportion), the distance is zero, whereas, the distance is one if completely different subpopulations.
Results Summary - The distance, the statistical significance (p-value), the subpopulations that are different in two profiles, and phylogenetic tree plot and tSNE plot to show the subpopulation structures.
Subpopulation view - heatmap of gene expression that mark the statistical significant subpopulations, and predicted cell types of these subpopulations based on Mouse Cell Atlas [Han et al., 2018, Cell, 1091-1107] if scRNA-seq are derived from mouse samples.
Subpopulation structure:
| Subpopulation | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
|---|---|---|---|---|---|---|---|---|---|---|
| S1 | 411 | 113 | 272 | 213 | 64 | 46 | 28 | 7 | 8 | 0 |
| S2 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 534 |
| Subpopulation | 10 | 1 | 3 | 4 | 2 | 5 | 6 | 7 |
|---|---|---|---|---|---|---|---|---|
| Proportion difference | 1.00 | 0.35 | 0.23 | 0.18 | 0.10 | 0.06 | 0.04 | 0.02 |
Figure 1: Phylogenetic tree plot of subpopulation structures. Cell subpopulations that are enriched in the first scRNA-seq profile labeled as red, while subpopulations enriched in the second profile are labeled as blue. The line width is proportion to the enrichment. Thicker lines suggest more enriched in one profile. Cell subpopulations that are similar in both samples are labeled as black.
Figure 2: tSNE plot of subpopulation structures. Each subpopulation that are different between two profiles are denoted as different colors, and subpopulations that are similar between two profiles are denotes as gray. Cells from the first profile are denoted by circle, while cell from the second profile are denoted by star.
Figure 3: Heatmap of gene expressions that mark the cell subpopulation 10
Figure 4: Predicted cell types of the cell subpopulation 10
Figure 5: Heatmap of gene expressions that mark the cell subpopulation 1
Figure 6: Predicted cell types of the cell subpopulation 1
Figure 7: Heatmap of gene expressions that mark the cell subpopulation 3
Figure 8: Predicted cell types of the cell subpopulation 3
Figure 9: Heatmap of gene expressions that mark the cell subpopulation 4
Figure 10: Predicted cell types of the cell subpopulation 4
Figure 11: Heatmap of gene expressions that mark the cell subpopulation 2
Figure 12: Predicted cell types of the cell subpopulation 2
Figure 13: Heatmap of gene expressions that mark the cell subpopulation 5
Figure 14: Predicted cell types of the cell subpopulation 5
Figure 15: Heatmap of gene expressions that mark the cell subpopulation 6
Figure 16: Predicted cell types of the cell subpopulation 6
Figure 17: Heatmap of gene expressions that mark the cell subpopulation 7
Figure 18: Predicted cell types of the cell subpopulation 7
devtools::session_info()
## setting value
## version R version 3.4.3 (2017-11-30)
## system x86_64, mingw32
## ui Rgui
## language (EN)
## collate English_United States.1252
## tz America/Chicago
## date 2018-05-16
##
## package * version date
## ape * 5.1 2018-04-04
## backports 1.1.2 2017-12-13
## base * 3.4.3 2017-12-06
## BiocInstaller * 1.28.0 2017-10-31
## bookdown 0.5 2017-08-20
## cluster 2.0.6 2017-03-10
## colorspace 1.3-2 2016-12-14
## compiler 3.4.3 2017-12-06
## curl 3.1 2017-12-12
## datasets * 3.4.3 2017-12-06
## devtools * 1.13.4 2017-11-09
## digest 0.6.15 2018-01-28
## evaluate 0.10.1 2017-06-24
## git2r 0.20.0 2017-12-17
## graphics * 3.4.3 2017-12-06
## grDevices * 3.4.3 2017-12-06
## grid 3.4.3 2017-12-06
## GUniFrac * 1.1 2018-02-12
## highr 0.6 2016-05-09
## hms 0.4.1 2018-01-24
## htmltools 0.3.6 2017-04-28
## httpuv 1.3.5 2017-07-04
## httr 1.3.1 2017-08-20
## kableExtra * 0.8.0 2018-04-05
## knitr * 1.20 2018-02-20
## lattice * 0.20-35 2017-03-25
## limma * 3.34.4 2017-12-14
## magrittr 1.5 2014-11-22
## MASS 7.3-47 2017-02-26
## Matrix 1.2-12 2017-11-20
## matrixStats * 0.52.2 2017-04-14
## memoise 1.1.0 2017-04-21
## methods * 3.4.3 2017-12-06
## mgcv 1.8-22 2017-09-24
## mime 0.5 2016-07-07
## miniUI 0.1.1 2016-01-15
## munsell 0.4.3 2016-02-13
## nlme 3.1-131 2017-02-06
## parallel 3.4.3 2017-12-06
## permute * 0.9-4 2016-09-09
## pillar 1.2.1 2018-02-27
## pkgconfig 2.0.1 2017-03-21
## plyr 1.8.4 2016-06-08
## questionr 0.6.2 2017-11-01
## R.methodsS3 * 1.7.1 2016-02-16
## R.oo * 1.21.0 2016-11-01
## R.utils * 2.6.0 2017-11-05
## R6 2.2.2 2017-06-17
## Rcpp 0.12.15 2018-01-20
## readr 1.1.1 2017-05-16
## rlang 0.2.0.9000 2018-03-07
## rmarkdown * 1.8 2017-11-17
## rmdformats * 0.3.3 2017-05-18
## rprojroot 1.3-1 2017-12-18
## rstudioapi 0.7 2017-09-07
## Rtsne * 0.13 2017-04-14
## rvest 0.3.2 2016-06-17
## scales 0.5.0.9000 2018-03-07
## scUnifrac * 0.9.4 2018-05-16
## shiny 1.0.5 2017-08-23
## splines 3.4.3 2017-12-06
## statmod * 1.4.30 2017-06-18
## stats * 3.4.3 2017-12-06
## stringi 1.1.6 2017-11-17
## stringr 1.2.0 2017-02-18
## tibble 1.4.2 2018-01-22
## tools 3.4.3 2017-12-06
## utils * 3.4.3 2017-12-06
## vegan * 2.4-5 2017-12-01
## viridisLite 0.3.0 2018-02-01
## withr 2.1.1.9000 2018-03-07
## xml2 1.2.0 2018-01-24
## xtable 1.8-2 2016-02-05
## yaml 2.1.16 2017-12-12
## source
## CRAN (R 3.4.4)
## CRAN (R 3.4.3)
## local
## Bioconductor
## CRAN (R 3.4.3)
## CRAN (R 3.4.3)
## CRAN (R 3.4.1)
## local
## CRAN (R 3.4.3)
## local
## CRAN (R 3.4.3)
## CRAN (R 3.4.3)
## CRAN (R 3.4.3)
## CRAN (R 3.4.3)
## local
## local
## local
## CRAN (R 3.4.4)
## CRAN (R 3.4.3)
## CRAN (R 3.4.3)
## CRAN (R 3.4.1)
## CRAN (R 3.4.1)
## CRAN (R 3.4.3)
## CRAN (R 3.4.4)
## CRAN (R 3.4.4)
## CRAN (R 3.4.3)
## Bioconductor
## CRAN (R 3.4.1)
## CRAN (R 3.4.3)
## CRAN (R 3.4.3)
## CRAN (R 3.4.1)
## CRAN (R 3.4.1)
## local
## CRAN (R 3.4.3)
## CRAN (R 3.4.0)
## CRAN (R 3.4.4)
## CRAN (R 3.4.1)
## CRAN (R 3.4.3)
## local
## CRAN (R 3.4.1)
## CRAN (R 3.4.3)
## CRAN (R 3.4.1)
## CRAN (R 3.4.1)
## CRAN (R 3.4.4)
## CRAN (R 3.4.1)
## CRAN (R 3.4.1)
## CRAN (R 3.4.3)
## CRAN (R 3.4.1)
## CRAN (R 3.4.3)
## CRAN (R 3.4.3)
## Github (tidyverse/rlang@68a84a6)
## CRAN (R 3.4.3)
## CRAN (R 3.4.4)
## CRAN (R 3.4.3)
## CRAN (R 3.4.3)
## CRAN (R 3.4.1)
## CRAN (R 3.4.3)
## Github (hadley/scales@d767915)
## Github (liuqivandy/scUnifrac@8a434d3)
## CRAN (R 3.4.1)
## local
## CRAN (R 3.4.2)
## local
## CRAN (R 3.4.2)
## CRAN (R 3.4.1)
## CRAN (R 3.4.3)
## local
## local
## CRAN (R 3.4.3)
## CRAN (R 3.4.3)
## Github (jimhester/withr@5d05571)
## CRAN (R 3.4.3)
## CRAN (R 3.4.1)
## CRAN (R 3.4.3)